Deconfinement and Dissipation in Quantum Hall "Josephson" Tunneling 
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The zero-bias tunneling resonance in quantum Hall bilayer systems is investigated via numerical 
simulations of the classical two dimensional XY model with a symmetry-breaking field. Disorder is 
included in the model, and is shown to nucleate strings of overturned spins proliferated through the 
system, with unpaired vortices and antivortices at their endpoints. This string glass state supports 
low energy excitations which lead to anomalously large dissipation in tunneling, as observed in 
experiment. The effect of an in-plane magnetic field is discussed. 
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Introduction. Bilayer quantum Hall systems have at- 
tracted much attention because, for small layer separa- 
tions and vanishing interlayer tunneling, they support 
states with spontaneous interlayer coherence. This bro- 
ken symmetry leads to a Goldstone mode [1] which may 
be described as a spin wave in an XY ferromagnet [2] . If 
one introduces interlayer tunneling, the system may al- 
ternatively be described as a large area Josephson junc- 
tion. In the XY language, the tunneling appears as a 
magnetic field tending to align the spins along the x 
direction- i.e. an in-plane field. While this symmetry- 
breaking field makes the physics of this system com- 
pletely different than the heavily studied model without 
an in-plane held, vortex deconfinement does occur, albeit 
not via the usual Kosterlitz-Thouless mechanism [3,4]. 
The Josephson junction analogy suggests these systems 
might support a Josephson effect [5]. Recently, a sharp 
conductance resonance near zero interlayer voltage bias 
was detected [6,7] which is highly reminiscent of the DC 
Josephson effect. The width of this resonance, however, 
is surprisingly large, and perturbative treatments [8-10] 
of the system have so far not been able to account for the 
dissipation. As we demonstrate below, the resolution of 
this mystery is likely related to the unbinding of vortices. 

In this article, we report on the results of extensive sim- 
ulations of this system, both with and without disorder 
(introduced via a random vector held described below), 
focusing on the effects of vortex deconfinement. Our re- 
sults have important implications for this model and its 
many realizations. Most significantly, when the disorder 
is sufficiently strong, a new type of low-temperature state 
emerges in which strings of overturned spins proliferate, 
with unpaired vortices at their endpoints. Low-energy 
excitations appear that are localized at the vortices, and 
the strings themselves behave in many respects like a 
pinned manifold. These observations are consistent with 
glassy behavior, so we call this a string glass state. The 
string glass is likely a "frozen-in" version of the decon- 
fined state found recently in studies of this system in 
the clean limit [3,4]. While we focus below on the im- 
plications of this state for the bilayer quantum Hall sys- 
tem, analogous states should also exist in planar Joseph- 
son junctions, real planar magnets, and two-dimensional 



crystals in an appropriate disorder environment, with 
consequences analogous to those we find for the quan- 
tum Hall system. 

For tunneling in the bilayer quantum Hall system, 
our results have a number of experimental implications. 
Firstly, the localized excitations near the vortices are 
directly excited by a tunneling current, leading to an 
anomolously large dissipation at low current. Such dissi- 
pation is known to exist in real systems [6], but its large 
magnitude so far has dehed explanation [8-11], partic- 
ularly its tendency to have a weak temperature depen- 
dence relative to what one might hnd if the dissipation is 
caused by quasiparticles. By contrast, we find that dis- 
sipation associated with the vortices in the string glass 
state intrinsically has a weak temperature dependence, as 
illustrated in Fig. 1. Above a critical current scale i c , we 
find that the strings depin and the resistance, as well as 
the voltage noise, is greatly enhanced. An observation of 
a peak in the noise near i c would give direct experimental 
evidence of the existence of the string glass state posited 
here. Finally, we show that an in-plane magnetic held Bu 
tends to align the strings along the held direction, and 
for large enough Bu the density of strings increases. At 
low temperature, the strings are pinned by disorder, so 
the zero-bias resonance remains in place. Again, this has 
been observed experimentally [7], but to our knowledge 
this is the first qualitative explanation of why this occurs. 

Simulation Model. To conduct our simulations we 
place our system on an L x L square lattice, whose lattice 
constant oq should be understood as the underlying mi- 
croscopic length scale of the system, the magnetic length 
£q = ^Jhc/eB, with B the perpendicular magnetic field. 
We impose periodic boundary conditions on this lattice. 
Our Hamiltonian is H = Hxy + Hp, with 
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= -K 



cos[0(r)-0(r')]-^cos0(r), (1) 
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A„> sin[0(r)-0(r')]. 



(2) 



In Eq. 1, 6(r) represents the angle of a planar spin lo- 
cated at lattice site r, J2<rr'> denotes a sum over near- 
est neighbors, and h represents the effect of a magnetic 
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field tending to align spins along 0(r) = 0, whose un- 
derlying source is the interlayer tunneling matrix ele- 
ment of the electrons. In Ho, A rr i is a random vector 
field residing on the bonds, which may be understood 
as follows. At long wavelengths, the state of the bilayer 
quantum Hall system may be specified by a three dimen- 
sional pseudospin field S(r). S z here represents the dif- 
ference in electron density between the two layers, and 
S x , S y are the components of the planar spin field. S z 
ideally vanishes when the electric potential of the two 
layers is perfectly balanced; in practice, there are im- 
purities which tend to push the electrons into one layer 
or the other. These same impurities will tend to locally 
lower or raise the total electron density, which can be 
accomplished by introducing a non-vanishing topologi- 
cal density q(r) = e^ u S ■ [<9 M S x <9„S]/87r, which in the 
quantum Hall context is proportional to the real charge 
density [12]. The topological density contains a term of 
the form A ■ V0 [9], with A proportional to the induced 
S z (r). It is this contribution to the disorder we model, 
taking d^d — ► sin[#(r) — 0(r')], with < rr' > the near- 
est neighbor bond along the (j, direction, to account for 
the finite grid spacing in the simulation [14]. We take 
A p (r) = ^2 r , e~' r ~ r ' ^^(r'), with drawn from a 
uniform distribution satisying —A < < A [15]. Be- 
cause disorder due to charged impurities in most quan- 
tum Hall systems varies slowly in space, the results we 
present are for relatively large values of £ p . Our qualita- 
tive results however also hold for smaller values of £ p . 

To simulate dynamics in this model, we adopt a 
Langevin equation of motion, 

d 2 9(r) 8H „ . d0(r) . 

r above represents an effective moment of inertia for an 
XY spin. Since 6 is conjugate to S z in the full quantum 
theory, T originates from the energy cost for creating 
charge imbalance across the layers, and is proportional 
to the capacitance per unit area of the bilayer. The £ 
and 7 terms are respectively a random torque and a uni- 
form damping, introduced to model coupling to a heat 
bath. To satisfy the fluctuation-dissipation theorem, the 
random torques are drawn from a distribution satisfying 
< CO, t)C(r',tf) >= 2jT5 r yS(t - t 1 ) with T the temper- 
ature of the system. Finally, the last term is a current 
which may be imposed to drive the system [13]. 

Results. We begin by studying the resistance of the 
system near zero bias. This can be conveniently accom- 
plished by examining the dynamics of the XY spins via 
the extended variables u(r, t) = J dt de ^ . The vari- 
ables u are identical to the original angles except they 
are unbounded. Their dynamics are analogous to those 
of particles in a washboard potential [16], and in the ab- 
sence of a drive current {id = 0) they diffuse. It is easily 
confirmed that < u 2 (r, t) >oc t for large enough time 
t, where < • > denotes an average over sites and disor- 



der configurations. We thus define a diffusion constant 
D = linii^ood < u 2 {r,t) > /dt. The "conductivity" of 
the w's is then given by the Einstein relation a = D/T, 
which quantifies their linear response to the drive cur- 
rent. Since uniform average motion of the it's generates 
a (Josephson) voltage, this means a is proportional to 
the tunneling resistance of the system, R. 
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FIG. 1. Tunneling resistance averaged over 16 disorder 
realizations and 10 6 Langevin sweeps for states in string 
glass state. Length, time and energy units chosen so that 
r = K = oo = 1; h^p — .25, 7 = 0.143 in these units. Inset: 
Corresponding resistance for no disorder. 

Fig. 1 illustrates some typical results for a measure- 
ment of this quantity in a series of runs. In these simula- 
tions, the system was annealed from above the Kosterlitz- 
Thouless temperature (Tkt — 0.892_ftT), with tempera- 
ture steps AT = 0.05. In the inset we show the result for 
a clean system, for which the tunneling resistance drops 
exponentially at low temperatures. At higher disorder 
levels (main figure) the low temperature resistance is in- 
creased by many orders of magnitude, demonstrating the 
crucial role played by disorder. 

To see why this is so, one may examine the detailed 
behavior of the XY spins. Fig. 2 illustrates a snapshot 
of a system configuration at low temperature. We show 
only spins for which cos6>(r) < 0; i.e., spins whose x 
component is opposed to the direction specified by the 
symmetry breaking field h. Qualitatively, one may see 
these form string-like structures. We find that narrow 
necks and endpoints of these strings are often associated 
with isolated vortices [3,17]. Also shown as circles are the 
100 most active sites, those with the greatest change in 
u(r) over the course of the simulation. These are clearly 
located in the vicinity of the isolated vortices, demon- 
strating that the barrier for overturning spins in such lo- 
cations is greatly reduced compared to other locations in 
the system. This is the underlying cause of the enhanced 
dissipation in these systems. At low temperatures and 
in the absence of a drive current id, we find that the 
strings (vortices) are pinned. Physically we expect this 
because these objects carry charge density (net charge) 
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and so may be held in place by disorder. Because of 
the quenched nature of the disordered state, we call this 
a string glass. Low-energy, localized excitations are a 
common feature of glasses [18], and their presence here 
supports our characterization of this as a glassy state. 
[19] 




FIG. 2. Snapshot from run with T = K = a = 1, L = 79, 
h^ 2 p = .25, 7 = 0.143, £ P /L = 0.13, A$ = 3.18. Only spins 
with cos 8 < shown. Circles indicate most active spins, 
which are localized near vortices (+) and antivortices (*) 

In the string glass state, a drive current acts as an ef- 
fective force on the strings. In a typical run, at very low 
drive the average "displacement" J d 2 r u(r,t)/L 2 for a 
single disorder realization has a small slope. As id is in- 
creased, one finds occasional short intervals for which this 
slope is increased, occuring when string sections become 
briefly depinned. With further increasing id these events 
become more frequent, until a critical drive current i c is 
reached at which the strings essentially move freely. In 
this situation, spins that were previously nearly static 
now rotate much more frequently as strings pass through 
them (inducing a Joscphson voltage), and the dissipation 
in the system is greatly enhanced [20]. This is apparent 
in Fig. 3, which illustrates the current- voltage character- 
istic for a single disorder configuration, with each data 
point averaged over eight different initial seeds for the 
random force. For each of these seeds we have averaged 
over 2 x 10 6 Langevin sweeps, with the first 25% of these 
thrown away for equilibration. The onset of "normal" 
dissipation in this system is apparently a dcpinning phe- 
nomenon. Because of the stick-slip motion in the vicinity 
of i c described above, there is considerable broadband 
noise in the dissipation for id ~ i c . This is typical of a 



dcpinning phenomenon, and should generate a peak in 
the noise spectrum of the power dissipated. An exper- 
imental observation of such a noise peak would confirm 
the string glass nature of the state, since it reflects the 
depinning of the strings. 

An in-plane component of the magnetic field leads to 
further interesting physics. For an appropriate choice 
of gauge, the in-plane field introduces a phase in the 
tunneling matrix clement which causes the XY spins to 
tumble spatially with a wavevector Q = {de/Tic)i x By, 
with d the separation between layers and By the in- 
plane component of the field. This can be incorporated 
into Hxy by making the substitution cos[#(r) — 6{v')\ — ► 
cos[6>(r) — 8{v') — Q Tr >\ in Eq. 1, where Q rr > = Qa resides 
on nearest neighbor bonds along the direction of Q. Note 
that because of the periodic boundary conditions, we are 
constrained to choose values of Q such that uq = LQ/2tt 
is integral. In the string glass state, for small values of 
uq, the main qualitative effect we observe in configura- 
tions produced by simulated annealing is a tendency for 
the strings to align along Q. For larger uq, we find the 
strings become narrower and the configurations appear 
more ordered, resembling a disordered soliton lattice [21] 
- with the strings being essentially identical to solitons. 

This qualitative transition is reflected in the critical 
current i Cl as illustrated in Fig. 3. One can see i c is 
little changed from its uq = value for uq = 1,2, 
whereas for uq = 3 and above, i c becomes noticeably 
suppressed. The existence of a parallel field scale above 
which the dcpinning current decreases can be understood 
in terms of collective pinning theory [22]. For large val- 
ues of Q, the soliton lattice may be thought of as an 
elastic medium, which is deformed and pinned by the 
disorder. This deformation involves a length scale L c for 
a typical domain size, determined by balancing the pin- 
ning and elastic distortion energies [22] . The length scale 
is meaningful provided L C Q > 2ir; if this is not satisfied 
then L c is smaller than the average distance between soli- 
tons and one enters a strong pinning regime. Following 
this reasoning [23] we obtain an estimate for the criti- 
cal wavevector above which one enters the weak pinning 
regime, Q c — hAy/n^/2nK 2 , where n p is the density of 
pinning centers. This implicitly defines a critical paral- 
lel field scale B c above which the depinning drive field i c 
will decrease with By . Below this, the strings are strongly 
pinned by disorder, and Bu has little effect. This is con- 
sistent with our simulation results, as illustrated in Fig. 
3. 

Our results show that the zero-bias tunneling reso- 
nance is not shifted to higher voltage by an in-plane field, 
consistent with experiment but in contrast to perturba- 
tive theoretical treatments [8,9]. The latter results sug- 
gest that a resonance should be seen at a voltage cor- 
responding to the Goldstone collective mode frequency 
huj(Q) = eV. A weak structure is observed obeying 
this relation [7]. While the simulations reported here 
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do not exhibit this behavior, it is likely that under differ- 
ent circumstances they would. Well above the depinning 
threshold, when sliding of the soliton lattice proceeds at 
a velocity such that a pinning center crosses a soliton 
once per period of a collective mode, the resulting corn- 
mensuration acts as an increased friction that lowers the 
tunneling resistance of the system. A perturbative anal- 
ysis [23] supports this, and suggests the width of such 
a resonance would be of the order AV ~ y/KVj/Te. 
For this to be observable in simulation, one needs very 
small values of 7, requiring impractically long equilibra- 
tion times. Such a small effective 7 is consistent with 
experiment. 
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FIG. 3. Induced voltage vs. id for different parallel mag- 
netic fields. Units chosen so that T = K = 00 = 1. 

Implications and Summary. The picture that emerges 
from these simulations, and the string glass state they 
suggest, has several experimental implications. These 
include the noise peak near the depinning transition and 
the critical parallel field discussed above, both of which 
can in principle be observed. In addition, one would ex- 
pect that small systems, for which the number of vortices 
generated by disorder is small, would generate strong 
sample to sample fluctuations in the tunneling resistance. 
Another signal would be the observation of a low tem- 
perature scale Tf « VKMo below which the isolated 
excitations freeze out and the tunneling resistance drops 
rapidly [23,24]. Tf is well below experimentally acces- 
sible temperatures for reported sample parameters [6], 
but might be observed in samples with large h and £0. 
Finally, since vortices carry charge e/2 [2] in this system 
and are deconfined, it is interesting to speculate that this 
fractional charge might be observed in shot noise exper- 
iments [25]. 
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